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Abstract 

The mass spectrum and the pionic decay widths of the P-wave D mesons are studied using the 
Dirac equation for the light quark. Our aim is to determine the Lorentz property of the confinement 
force. We find that the Lorentz scalar confinement is consistent with the mass spectrum, while a 
significant mixture of the Lorentz vector confinement is necessary to explain the decay widths. 
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I. INTRODUCTION 



Confinement of color is the most outstanding problem in hadron physics. Quantum 
chromodynamics (QCD) has revealed that the ground state of pure gluonic QCD is in the 
confined phase, and that static (heavy) quarks are confined as is shown by the area law 
behavior of Wilson loops. ^| Yet, confinement of dynamical quarks is not well established, 
although we believe that light quarks are confined as well and all hadrons are white. 

In phenomenology, on the other hand, various well-tuned quark models are available in 
explaining hadron spectra. The simplest system may be the heavy quarkonia, such as cc 
(J/i/j, ip' ■ ■ ■ ) and bb (T). Their spectra are so simple as to be reproduced by a nonrelativistic 
Hamiltonian quark model. ^ The Cornell potential, 3] which consists of a linear confinement 
potential and a short-range color-Coulomb force, reproduces overall level structure and also 
agrees very well with the lattice QCD static potential. The recoil (or relativistic) corrections, 
up to the order 1/Mq, can be computed for the single gluon exchange process and are shown 
to give reasonable spin dependences for fine structures of the spectra. j^J 

One important question, which has not been properly answered in the lattice QCD, 
is the covariance property of the confinement. In the Hamiltonian quark model, as an 
approximation of the full QCD dynamics, the inter-quark potential may have various Dirac 
structures at the vertex. For instance, the one-gluon exchange potential is assumed to have 
a 7 M at each vertex. In the static potential limit only the time component 7 survives, which 
of course give a Coulomb potential. On the other hand, if we exchange a few gluons between 
the quarks, for instance, the Dirac structure at the vertex can be further complicated. 

In addressing the above problem, we choose a semi-simple system, that is, bound states 
of a heavy antiquark (Q) and a light quark (q). Concretely, we treat here the D- meson 
spectrum, which contains a charm (anti) quark and a light quark. This heavy-light system 
has a simple limit if the heavy quark mass Mq is infinite. The system reduces to a single 
quark moving under a potential which is independent from the heavy quark spin. The spin 
independence is a result of the heavy quark symmetry of QCD. 5] Now the light quark is 
supposed to move relativistically and therefore we need to treat this system (in the large 
Mq limit) in terms of the Dirac equation. [f| 

Such Dirac equation is required to have a confinement potential, but the Lorentz property 
of the confinement potential is not a priori given. Without derivative (or velocity depen- 
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dence), we have two choices in general. Either a Lorentz scalar potential, which appears in 
the Dirac Hamiltonian as /3S(r), or the time component of a Lorentz vector potential, given 
as V(r). Namely, the Dirac Hamiltonian of the form 



H = a ■ v + + /3S(r) + V(r) (1) 

is supposed to represent the system. Besides, we also need the color-Coulomb part of the 
interaction for this system as well, which is attributed to the one-gluon exchange interaction 
and is therefore included in the Lorentz vector part V(r). 

The purpose of the present paper is to study how we can determine the Lorentz property 
of the confinement by using the Qq meson observables, i.e., the masses and decay amplitudes. 
For this purpose we employ a linear confinement potential whose total strength is fixed but 
the ratio of the Lorentz scalar and vector components is allowed to vary, 

V coni = f3(l - v)br + vbr (2) 

The parameter v denotes how much of the confinement is regarded as the Lorentz vector. 
By varying v, we study the D meson spectrum and decays and attempt to determine which 
v is the most appropriate. We are particularly interested in the P-wave spectrum of the D 
mesons, because there the spin-orbit interaction plays an essential role. It is well known that 
the Coulombic potential in the Dirac equation gives a splitting of the pi/ 2 and p 3 / 2 states, 
where the p%/2 is lower in energy than ^3/2- The order of the two p states is determined by the 
sign of the effective spin-orbit potential, and depends on whether the potential is regarded 
as the Lorentz vector or scalar. [7] Indeed this splitting is reversed when the Dirac equation 
for a pure scalar attractive potential is solved regardless of the shape of the potential. The 
MIT bag modelQ] is an example, where the ground lsi/2 is followed by the 1st excited state 
2p 3 / 2 instead of 2pi/%- It is therefore extremely interesting to find the order of the 2pi/ 2 and 
2j?3/2 states in the D meson spectrum. 

In reality, however, as the heavy quark spin is not completely decoupled, we need to 
take into account 1/Mq correction terms and especially the hyperfine splitting due to the 
spin-spin interaction. 

In Sect, m we present the Dirac equation and its solutions relevant to the discussion in 
this paper. In Sect. Hill 1/Mq correction terms are introduced both for the scalar potential 
as well as the vector potential. In Sect. IIVI we calculate the pionic decay widths of the D 



mesons. The decay widths are computed in the first-order perturbation theory and agree 
with the results from the heavy quark symmetry. In Sect. 13 we first explain how we 
determine the values of the parameters and then present the results of our calculation. In 
Sect. IV11 a conclusion and discussion are given. 



II. WAVE FUNCTIONS 



We start with the simple picture of the D mesons, i.e., a light quark bound in a potential 
created by a heavy anti-quark, which sits still at the origin. We employ the Dirac equation 
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with a linear plus Coulomb potential to describe the motion of the light quark. 

As shown in Eq. (J2J), we assume that the linear potential may include scalar and vector 
components. 



C = 



i$ — m + (3 br{(l — v) + (3v\ 

r 



* (3) 



where m is the constituent mass of light quark, a = ~a s , and b is the string tension. The 
corresponding Hamiltonian is given by 

H = a. ■ p + (3m — — + (3(1 — v)br + vbr (4) 

The total angular momentum j is conserved for the central potential in Eq. (jljl. and then 
the angular part of the wave function is written by the spinor spherical harmonics, 

yjm(M) = J2 (£m e sm s \jm}Y £me (6,(f))x ms (5) 

mi,m s 

We parametrize the wave function as 

where I and £' denote the orbital angular momentum of the large and small components 
respectively. Then the equations of motion for the radial parts are given by @ 

a „\ „, . I d k 



(m + (l- v)br + vbr - - - e) F(r) = I- ) G(r) 

\ r ) \dr r J 

(-m - (1 - v)br + vbr-^~ E^j G(r) = (-^- - ^ J F(r) (7) 



where k = ±(j + |) for j = ±|. Explicitly, k — — 1 for Isi, k = —2 for 2p3 and k = 1 for 
2pi. E is the energy eigenvalue, for which we only consider the positive energy solutions. 



Wave functions and energy eigenvalues are acquired by solving this equation. Some of the 
obtained wave functions are shown in Fig. ^ We will use them to calculate the 1/Mq 
correction to the mass spectrum as well as the decay widths in the following sections. 
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FIG. 1: Wave functions of ls±, 2ps and 2pi states. (m = 300MeV |a s = 0.360 b = 

222 6 

0.190GeV 2 v = 0%) 

III. 1/Mq correction terms for the d meson masses 

Once we solve the Dirac equation for the light quark in the linear plus Coulomb potential, 
we are ready to compute the mass spectrum of the D mesons. In the heavy quark limit, 
the heavy quark spin is decoupled and we have sets of degenerate states. Among them, 
here we are interested in D (J n = 0~) and D* (J n = 1~) states from the lsi eigenstate 
of the light quark, D* (J n = 0+) and D{ (J* = 1+) from 2pi, and D 1 (J n = 1+) and D* 2 
(J 71 " = 2 + ) from 2j?3. In the real world, the charm quark is not infinitely heavy so that we 
have splittings in each pair. The most popular one is the spin-spin splitting in the S states, 
or D — D* splitting. The splittings are commonly attributed to the interaction of higher 
orders in 1/Mq expansion, where Mq is the heavy quark mass. It is also important to note 
that the Di and D\ states may mix with each other when the heavy quark spin is allowed 
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to flip. Thus, the 1/Mq higher order terms will mix D\ and D\. 

We here derive the 1/Mq interaction terms associated with the linear plus Coulomb 
potential, which are regarded as relativistic corrections coming from the Lorentz scalar 
and/or vector couplings. The corrections for the Coulomb part is similar to the Fermi-Breit 
interaction for the heavy Q — Q system. We consider the one-gluon exchange diagram in 
the Coulomb gauge, which gives the gluon propagator 

47T 

and compute the qQ scattering (invariant) amplitude to the leading order, 
M = -^(u[Y Ul )D, u (q)(u' 2 Yu 2 ) 

= ~ {K 7 V)K7%)A,o + («iy«i)(tiW*tt2)A*} (9) 

Here, u± and u 2 denote the plane wave Dirac spinors of the light and heavy quarks, respec- 
tively. Then, we replace the heavy quark spinors in Q by their static expansion, 



, 2(P H * J*l* x | do) 

<2E Q \ JE Q -M Q n-a X 



and the amplitude up to the order to 1/Mq is given by 

M = -ajtwaW) ^00 + -^-(u'^Uj)^ ( itr x (-q) + 2p 2 - q^Afc} (H) 

1st part 2nd part 3rd part 

Now, we determine the effective potential (which is in general, velocity dependent) so 
as to recover this amplitude in the Born approximation. The potential corresponding to 
the first part is the Coulomb potential, —a/r, which is of order Mq. Energy eigenvalue 
corresponding to this is already acquired by solving (JJJ). The second and third parts induce 
1/Mq terms. The potential corresponding to the second part is given by 



Mi 



a 



( So x 4 



V S Q x ^ 



(12) 



Because this term may change the spin of the heavy quark, this leads to a mixing between 
Di and D*. The potential corresponding to the third part is given by 



a 1 
Mq27 



1 <t ■ p H — r • p 



cr ■ r 



(13) 



Similarly, the 1/Mq terms generated by the linear potential can be obtained. There we 
need to distinguish the Lorentz scalar part and the Lorentz vector part of the linear potential 
and calculate the contribution of each term separately. 
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Finally, one more 1/Mq correction comes from the kinetic energy (recoil) of the heavy 
quark, i.e., 

p 2 /2M Q (15) 

where p is the momentum of the heavy quark. 

The 1/Mq Hamiltonians ()12)) . (jlH|) . (jl4j) and (|15jl are treated in the first order perturbation 
theory. As a result we obtain the masses of D* and D for the ground states, and those of D2, 
Di, D\ and Dq for the first excited states. Here, we consider the mixing between two states 
Di and D3, and obtain D\ and D\ as physical states. The detailed formulas are given in 
Appendix. |XJ 



IV. PIONIC DECAY WIDTHS OF D MESONS 

The pionic decay widths of the D mesons are calculated by using PCAC relation. [9] The 
decay widths into Dtt states are given 10] by 

T(D n ^Dn)= 3 , |q ^ \X DDn \ 2 (/ w = 93MeV) (16) 

where q is the momentum of the emitted pion, and the amplitude X^Dn ^ s given by 

X DDn =Jdr ¥ D (r) U - A • ^„(r)e-^ r (17) 

Here we neglect the final state interaction. Taking the direction of q along the z-axis, we 
obtain 

00 . 

X DDn = ^(-z)y(2£+l)(47r) / dr v&L(r)(75 - ^ Dn (r)j e ( q r)Y e °(r) (18) 

From the conservation of parity and angular momentum, one can easily see that D*, Dq, and 
D\ can decay into Dn, while D\ and D\ cannot. The explicit forms of XDD n {D n = D*,Dq 
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and D 2 ) are given by 
X DD *{q) 

X DD *{q) 



-|F (r)| 2 + l|G (r)| 2 } jo(gr) + ^|G (r)| 2 j 2 (gr) 



dr 



ziF^G^-G^F^jUqr) 

-t{F (r)F 1 (r)+G (r)G 1 (r)}j 2 (qr)\dr 
X DD * 2 (q) = J [V2i{F (r)F 3 (r) - ±G (r)G 3 (r)} ji(gr) 



W2i{F (r)G 3 (r) - G (r)F 3 (r)} j 2 (qr) - |-G (r)G 3 (r) |> j 3 (qr) 
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(19) 



fir 



The radial wave functions for the lsi, 2pi and 2j?3 states are labeled by 0, 1 and 3, respec- 
tively For the decays into the D*n final channel, we use helicity amplitude X^i*^ n , where 
the ho* and I%D n denotes the z-component of the spin of initial D n (final D*) meson. Since 

the z-component of the angular momentum is conserved, we obtain 

3 



D*7l) 



1^1 f I y++ |2,|y00 12,1^ |2\ 

■ , ^ \IAd*aJ +I^d*d„I +I^d*d„I / 



(20) 



87r/ 2 2j n + 

Now, Z)^, -Di and D\ can decay into .D*7r, while Dq cannot. Using the matrix elements given 
in Eq. (fT3j) . the decay width of T(D* 2 — > D*n) is given by 

3 3 



T(D 2 



D*n) 



2 8tt// 5 



(21) 



We consider the mixing of D\ and -D x according to Eq. (|A5|) . Denoting the eigenvectors of 
Eq. (|A5|) by (a, — /3) and (/3, a), we obtain 



r(£>r^L>*7r)(g) 



The heavy quark symmetry 0| predicts the relations 



3|q| 3 

3|q| 3 
8^ 



-^=aX DD »(g) - (3X DD *(q) 
-^=[3X DD *(q) + aX DD *(q) 



(22) 



5 v 5 



-rm 

5 v 



D*7l) 



(23) 



T(D* 2 - Dtt) 

r(z>* -> d*it) 

T(D* -> Dtt) = T(Zh 

One easily sees from Eqs. (|21j) and (|22|) that these relations are satisfied if we use the same 
|q| for all the decay widths. In reality, due to the phase space differences, these relations 
are not exactly met. 

When we later compare calculations with experiment, we use the observed total decay 
widths of the D mesons, because the pionic decay is dominant in the D meson decays. 
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V. DETERMINATION OF PARAMETERS AND RESULTS 



The Dirac equation has five free parameters; the constituent light quark mass m, the 
strong coupling constant a s , the string tension b, the constituent heavy quark mass Mq, and 
the zero point of energy m c °. The ambiguity of the zero point corresponds to uncertainty of 
m c ° in (EH), (H2J) and (IaH) . 

These parameters are to be determined phenomenologically with some conditions im- 
posed. First, the string tension is predicted by the lattice calculation so as to be 
0.17 ~ 0.20GeV 2 . The constituent charm quark mass should be about the half of the 
mass of J/ty, which is about 1.5GeV. We will leave out unphysical region of parameters, 
which are too different from these values. 

Now, these parameters are determined by using the well-determined masses of the D 
meson family, which are D, D*, Di and D^. Paying attention to the differences of the 
masses of the D mesons instead of their absolute values removes ambiguity of the zero 
point. Then, four parameters remain. But, the differences of the masses of D mesons give 
only three constraints. So, we cannot fit all the parameters. Therefore, we choose the 
constituent light quark mass m as 300MeV. The difference between D and D*, which is 
141.45MeV, determines Mq through (^— + bv)) in (jAl|) . After that, we determine the 
values of a and b, by taking into account the experimental values of D\ — D* mass splitting 
(413.85MeV) and of D* 2 - D x mass splitting (36.8MeV) . 

The parameter dependences of various quantities are shown in Figs. |2K From Fig. 
one sees that in order to reproduce the Di — D* mass splitting we must use a large value 
of b about 0.26GeV 2 . However, this value is not acceptable from the lattice calculation. 
Furthermore, in the region of ^a s > 0.35 by taking such a large value of b, we are forced 
to use too large value of Mq c in order to reproduce the experimental value of D* — D mass 
splitting. (See Fig. 0) 

Thus we give up reproducing the D\ — D* mass splitting, and employ b = 0.19 GeV 2 , 
which is within the lattice prediction. In this case the predicted D\ — D* mass splitting is 
about 340MeV, which is 70MeV smaller than the experimental value. 

From Fig. |21 one notices that small (large) a s is favored for small (large) v. However, 
because of the constraint on the value of Mq c shown in Figs. 0] and [7[ we cannot let a s too 
small (large). We therefore choose a = |a s within the region 0.32 < a{= |a s ) < 0.40 in 
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FIG. 2: Difference between the experimental data of D\ — D\ and our calculation for various a s 
and v. 
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FIG. 3: Difference between the experimental data of D\ — D* and our calculation for various a s 
and v. 
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FIG. 4: The constituent charm quark mass for various a s and v. 
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FIG. 5: Difference between the experimental data of D\ — D\ and our calculation for various a s 
and b. 
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FIG. 8: The masses of D mesons in P-wave. We fit m(-D|) = 2459MeV. The value of parameters 
are |a s = 0.36, b = 0.19GeV 2 and m = 300MeV. 
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FIG. 9: The mixing angle of D\ and D*. The value of parameters are |a s = 0.32, b = 0.19GeV 2 
and m = 300MeV. 



order to reproduce the experimental D* — D mass-splitting. The value of m c ° is fixed by 
m{D^) = 2459MeV, which is the most reliable among the P-wave states. 

The results are shown in Fig. |H| This figure indicates that D\ does not cross Di in the 
case of 0% < v < 60%. From Fig. |Hl one sees that the mixing angle of D\ and D\ decreases 
as v increases, while the difference between the eigenvalues of M{D3_) and M(Di) decreases 
as v increases. As the mixing angle is small, D\ is dominated by Di, and D\ is dominated 
hy Di in all region of v. This is caused by the off diagonal elements of the mixing matrix 
given by (jA4j) . As is seen from Fig. 1, the F 3 (r)Gi(r) product in Eq. (jA4|) changes its sign 
at around r ~ 3.0 GeV -1 . Because the Coulomb term is dominant at short distances and 
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FIG. 10: The pionic decay widths of D\ and D\. The value of parameters are |a s 
0.19GeV 2 and m = 300MeV. 
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FIG. 11: The pionic decay width of D\. The value of parameters are |a s = 0.32, b = 0.19GeV 2 
and m = 300MeV. 
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FIG. 12: The pionic decay width of Dq. The value of parameters are |a s = 0.32, b = 0.19GeV 2 
and m = 300MeV. 
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the linear potential is dominant at large distances, they tend to cancel when v increases. 

The pionic decay widths of D\ and D\ are more sensitive to the mixing effect than their 
masses, because the difference between T(Dz — > D*) and T(Di — > £)*) is very large. The 
Di(Ds) decaying into D*tt occurs mainly in the S-wave (D-wave), so the decay width is large 
(small). Then the difference of the decay widths T(Di — > D*) and Y(D\ — > D*) decreases 
as the mixing angle is increased. 



We compare our results with experimental data from Belle [11 1, CLEO|l2|], and Particle 
Data Groupj.ljj. M{D*) is larger than M{D{) in the CLEO data, while the order is reversed 
in the Belle data. Our results agree with the former. In the case of small v, the masses of 
the P-wave D mesons are consistent with the data of CLEO. But, the difference of the decay 
widths T(Di — > D*) and T(Dl — > D*) is smaller than the observed value. This problem 
is settled by increasing v according to the above argument. Choosing v = 40%, we can 
reproduce T(Dl — > D*)/T(D 1 — > D) rather well. But, in this case, P-wave spectrum does 
not reproduce the experimental data. Furthermore, the charm quark mass is required to be 
heavy, ~ 2GeV. 

The results are summarized as follows. 

1. The mass spectrum of the P-wave D mesons indicates that the confinement force is 
almost scalar or, even if vector component exists, it is about 10 percent at most. 

2. On the other hand, the pionic decay widths indicate that the confinement force includes 
significant amount of vector component. 



VI. CONCLUSION AND DISCUSSION 



The heavy-light quarkonium spectra are suitable for the study of the Lorentz property 
of the color confinement force. Recent new data of the P wave spectrum of the D meson 
systems are quite important and powerful in determining the spin-dependent interactions 
among quarks. We concentrate on the difference of the spin-orbit splitting from the Lorentz 
scalar and vector quark confinement, and attempt to determine how much Lorentz vector 
confinement is necessary to account for the P wave spectrum of the D mesons. 

For this purpose, we start with the heavy quark limit in which the motion of the light 
quark is treated by a single particle Dirac equation. The light quark is bound by a potential 
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that consists of the confinement and one-gluon-exchange pieces. By varying the vector 
component of the confinement, we study the spectrum as well as the pionic decay widths of 
the D mesons. 

We pick up three choices of v ; A: v — (pure scalar confinement), B: v = 0.1, and C: 
v = 0.4. The results are summarized in Table. HJ A problem common to all the parameter 
choices is that the splitting between IS* and 2P states is too small. We have found that 
unless we change some unnaturally large confinement potential, we are not able to reproduce 
the difference. Thus we give up fitting the parameters for this splitting. 

When we tune up the parameters mostly to the P-wave spectrum, then we obtain set A 
and set B, while the decay widths happen to prefer set C. The set A and B are consistent 
with the mass spectrum obtained at CLEO, but the pionic decay widths are not satisfactory, 
though the set B is slightly better. We note, however, that the recent data from Belle group 
predict the D\ mass much smaller than that given by CLEO, and in fact the order of D\ 
and Di is reversed. Within our searches of the parameters, we never see the reversal of this 
J = 1 states. It is highly desirable to determine the mass of D* precisely in order to reduce 
the vector amount of the confinement. 

The set C is unique as it predicts significant amount of Lorentz vector component in the 
confinement. The results are consistent with the pionic decay widths, but show a worse fit to 
the mass spectrum, and also the charm quark mass becomes large. As the experimental data 
for the decay widths of D\ and D\ are in good quality, it seems preferable to use these data 
to fix the vector amount of the confinement. We then see that a larger vector confinement 
gives a better fit. At the same time we find that the mixing of pi/2 and P3/2 states in the 
J = 1 mesons, Di and D\, decreases when the vector part of the confinement increases. This 
causes the pionic decay width of D\ to get lower, as the pure pz/2 state cannot decay into 
7r and S state by keeping the conserving the angular momentum and parity simultaneously. 
Thus the experimental data of the pionic decays suggest that the mixing is suppressed. In 
our approach this leads to the conclusion that the vector component of the confinement 
force must be significant in order to reduce the mixing strength. 

Similar approaches have been taken by Ebert et a/.Q] and also by Di Pierro et q/.jisl. 
The former considered the Qq systems in terms of the Dirac equation with 1/Mq correc- 
tions. They considered mixing of Lorentz vector confinement as well as the anomalous 
color-magnetic coupling. They fixed the ratio of the Lorentz vector and scalar confinement 
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Our results (m = 300MeV b = 0.190GeV 2 ) 


Experimental data 




A 
v = 
^ = 0.32 


B 

v = 0.10 
^ = 0.34 


C 

v = 0.40 
= 0.40 


Belle|llj 


CLEO[12l 


PDG|13J 


M(D*) [MeV] 


2459 


2459 


2459 


2460.7 ±2.1 ±3.1 




2459 


M{D*) [MeV] 


2459 


2458 


2435 


2400 ± 30 ± 20 


246ltg 




M(Di) [MeV] 


2426 


2421 


2411 


2423.9 ± 1.7 ±0.2 


2422.2 ± 1.8 


2422 


M(D*) [MeV] 


2356 


2357 


2334 


2290 ± 22 ± 20 






M(D*) - M(D t ) - 36.8 [MeV] 


-4.289 


1.165 


11.36 








Af (Di) — M(D*) — 413.85 [MeV] 


-80.46 


-97.83 


-105.3 








M Qc [MeV] 


1095 


1351 


1979 








r^-Dffl.^ [MeV] 


67.97 


65.89 


57.36 


46.4 ±4.4 ±3.1 




25t 8 7 


Fd* -D*7r [MeV] 


233.9 


276.2 


319.5 


380 ± 100 ± 100 


290t£° 




r Dl - D « w [MeV] 


75.30 


51.67 


28.06 


26.7 ± 3.1 ±2.2 




26 ±8 


r D ._ DT [MeV] 


360.8 


384.6 


415.0 


305 ± 30 ± 25 






r D ,_ D7r [keV] 


194.9 


187.1 


158.1 






142.12 ±0.07 


Mixing angle of D 3 and D 1 
[degree] 


13.64 


8.455 


3.178 









TABLE I: Our results and the experimental values of the masses and the pionic decay widths of 



the D mesons. 

as 2 : — 1. Their parameters are totally different from ours and their spectrum has a very 
high J = 1 state, which disagrees with current experimental data. 

Ref. |l5| studied the Dirac equation for the Qq system, but considered only the scalar 
confinement. Their results again show a high J = 1 state, while the Is — 2p splitting is 
reasonable. But their a s and b seem too large. 

In conclusion, it has been clarified how the masses and decay widths of the D mesons 
behave under the change of the Lorentz properties of the confinement force. It is interesting 
that the current experimental data of the decay widths prefer a large mixing of the Lorentz 
vector confinement. We, however, have some difficulties in fitting all the spectrum. We 
stress that the determination of the D\ mass and width is very important to fix the picture 
completely. Future experimental and theoretical studies are desirable. 

APPENDIX A: FORMULAS OF D MESON MASSES 

We present the explicit forms of the D meson masses up to O(jj^). Derivatives from the 
momentum operator are rewritten by means of the Dirac equation ((7|). From the lsi state 



16 



of the light quark, we have two states: 



|£>) = |J = 0(lai)> |£>*> = |J=l(l a i)> 



Their masses are given by 



m(D* 
m(D) 



m c ° + E + 

7T~ dr 
Mn J 



drF (r)G (r) \^ + bv 



/_2 N 

™ + &) {F (r)} 2 - (™ + 6(1 - 2v) ) {G (r)} 2 



bv 



a T \dr i(^ + i?o)({Fo(r)} 2 + {Go(r)} 2 ) + 



F (r)G (r) 



bv 



(mr + 6r 3 ) {F (r)} 2 - (mr + 6(1 - 2w)r 3 ) {G {r)} 
Jdrr(l + £ ) (Vo(r)} 2 + {G (r)} 2 ) 
l —Jdr {^m + (1 - t;)6r} 2 - jt;6r - °- - E }^ ({F (r)} 2 + {G (r)} 
l —jdr 2(1 -v)bF (r)G (r) (Al) 



2M 



From the 2j>3 state of the light quark, we have two states: 



|£>5> = |J = 2(2pa)> |£>s> = |J=l(2ps)> 



Their masses are given by 



m(D* 2 
m(D3 



m c ° + E 3 + 



dr 



Mr 



( 



4 
V 3 

m 



drF 3 (r)G 3 (r) + bv) 



+ 6 {F 3 (r)} 2 - 



m 



b(l-2v)){G 3 (r)} 



M c 
bv 



dr 



r 



— I 

M Q J 

'"' J dr [(mr + 6r 3 ) {F 3 (r)} 2 - (mr + 6(1 - 2w)r 3 ) {G 3 (r)Y 
Jdrr^ + E^ ({F 3 (r)} 2 + {G 3 (r)} 2 ) 



M, 



Q 



l^Jdr (^m + (1 - v)br} 2 - {vbr - °- - £ 3 }^ (V 3 (r)} 2 + {G 3 (r)} 
l — j 'dr 2(1 -v)bF 3 (r)G 3 (r) (A2) 
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From the 2pi state of the light quark, we have two states: 



| J Di) = |J = l(2pi)) |A)) = |J = 0(2pi)) 



Their masses are given by 
m(Di) 

2 

m(D*) 



m c ° + E 1 + 



-2 
(m 



M, 



Q 



drFAr^GAr) 



+ bv 



M, 

a 



rl dr [{v +b ) {Fl(r)}2 - (7 + 6(1 - 2v) ) {Gl(r)}2 



Mq 
bv 

Mq 
bv 



dr 
dr 



F 1 (r)G 1 (r) 



M, 



Q 



(mr + 6r 3 ) {^(r)} 2 - (mr + 6(1 - 2w)r 3 ) {Gi(r)} : 
Jdrr^ + E^j ^{ Fl (r)} 2 + {G l (r)Y 



■-^ I dr (^m + (1 - t;)6r} 2 - j^r - ^ - ^ (V\(r)} 2 + {(^(r)} 2 
■^L Jdr2{l-v)bF l {r)G l {r) 



(A3) 

The radial wave functions for the lsi, 2 01 and 2»3 states are labeled by 0, 1 and 3, respec- 

22 2 

tively. The off-diagonal element of the mixing matrix of Di and D3 is given by 

V2 1 



3 M, 



-/ rfr (F 3 (r)Gi(r) + F 1 (r)G 3 (r)) + bv^j (A4) 



Then the eigenvalues are determined by the secular equation given by 

m(Di) - A 

m mix m(Ds) - A 
Solving this equation, the masses of D\ and Di are given by 







(A5) 



m(Dl) 



m(Di) 



m(Di) + m{Ds) + \j (m(Di) — m(Ds)^j + m mix 



m 



Di) + m(D^) — J (m(Di) — m(Ds_f) + m D 



(A6) 



It is conventional to assign D* to the state with larger Di component. In our case, D{ 
happens always to be heavier than D ± . 



[1] G. S. Bali, Phys. Rev. 343 (2001), 1. 



18 



[2] S. N. Mukherjee et al, Phys. Rev. 231 (1993), 201. 

[3] E. Eichten et al, Phys. Rev. Lett. 34 (1975), 369.; Phys. Rev. D17 (1978), 3090; D21 (1980), 
203. 

[4] A. De Rujula, H. Georgi and S. L. Glashow, Phys. Rev. D12 (1975), 147; N. Isgur and 
G. Karl, Phys. Lett. B72 (1977) 109; Phys. Rev. D18 (1978), 4187; D19 (1979), 2653; D20 
(1979), 1191. 

[5] N. Isgur and M. B. Wise, Phys. Rev. Lett. 66 (1991), 1130. 

[6] M. Kaburagi et al, Phys. Lett. B97 (1980), 143; Z. Phys. C9 (1981), 213; J. Morishita, 

M. Kawaguchi, and T. Morii, Phys. Rev. D37 (1988), 159. 
[7] N. Isgur, Phys. Rev. D57 (1998), 4041. 

[8] A. Chodos et al, Phys. Rev. D9 (1974), 3471; D10 (1974), 2599. 

[9] S. Weinberg, Phys. Rev. 177 (1969), 2604. 
[10] D. Becirevic and A. Le Yaouanc, J. High Energy Phys. 03 (1999), 021. 
[11] Belle Collaboration, K. Abe, et al, ICHEP02 abstract 724, BELLE-CONF-0235. 
[12] CLEO Collaboration, S. Anderson et al, Nucl. Phys. A663 (2000), 647. 
[13] Particle Data group, Phys. Rev. D66 (2002), 010001. 
[14] D. Ebert, O. Galkin and R. N. Faustov, Phys. Rev. D57 (1998), 5663. 
[15] M. Di Pierro and E. Eichten, Phys. Rev. D64 (2001), 114004. 



19 



